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Abstract 



A condensed survey of various aspects of LEP 1 physics, e + e~ — > ff, is presented, with 
a critical examination of the ingredients that are available for theoretical predictions. A 
prototype of comparisons for 2f calculations can be found in hcp-ph/9902452| . [] 



1 Prolegomena 

Consider the LEP 1 measurements of c^had and A 1 as an example to discuss the data analysis strategy: 
each LEP experiment extracts a set of 

po d ^{M z) T z ,R e ^ T) A F ^}, (i) 

from their measured a and ^4 pB (realistic observables) . The four sets of PO are combined, taking correlated 
errors between the LEP experiments into account, in order to obtain a LEP-average set of PO. The latter 
is then interpreted, for example within the frame-work of the Minimal Standard Model. The practical 
attitude of the experiments is to stay with a Model-Independent (MI) fit, i.e. 

from RO —> PO a SM remnant (2) 

for each experiment, and these sets of PO are averaged. The extraction of Lagrangian parameters, 

M z ,m t ,M H ,a s (M 2 z ), and a(M 2 z ), (3) 

is based on the LEP-averaged PO. Since PO are determined by fitting RO, one has to clarify what is 
actually taken from the SM, such as imaginary parts, parts which have been moved to interference terms 
and photon-exchange terms making the MI results dependent on the SM. 

An alternative way to extract SM parameters from the LEP measurements is to average directly the 
SM parameters which have been obtained by each experiment through a SM fit to its own RO. This 
alternative should also be pursued by the experiments. 

The 2f emergency kit is illustrated in Fig. [j]. 



2 Exegesis of comparison 

All experimental informations are passed through decoders like TOPAZO and ZFITTER. Therefore, it is of the 
upmost importance to critically upgrade and to compare the complete TDPAZO and ZFITTER predictions || 
with a particular emphasis on demanding the following criteria: 
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1 The numerical results are based on the work of TOPAZO |j| and ZFITTER Q teams. At present, the following physicists 
are active members of the two teams: TOPAZO: G. Montagna, O. Nicrosini, G. Passarino and F. Piccinini; ZFITTER: D. 
Bardin, P. Christova, M. Jack, L. Kalinovskaya, A. Olshevski, S. Riemann, T. Riemann. 
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• Comparisons, after the upgrading, should be consistently better than what they where in earlier 
studies. 

• At the peak all relative deviations among total cross-sections and absolute deviations among asym- 
metries should be below 0.1 per-mill. 

• At the wings, typically yfs = M z ± 1.8 GeV, they should be below 0.3 per-mill. 

It is important to observe that today's comparisons are at the level of 10~ 4 , never been attempted 
before. 

3 PO, RO and all that 

Within the context of the SM the RO are described in terms of some set of amplitudes 

A SM — A 7 + A z + non-factorizable, (4) 

where the last term is due to all those contributions that do not factorize into the Born-like amplitude. 
Once the matrix element A SM is computed, squared and integrated, a convolution with initial- and final- 
state QED and final-state QCD radiation follows: 

<r(s) = J dz H in (z,s)H Rn (z,s)a(zs), (5) 

After de-convoluting RO of QED and QCD radiation aN set of approximations transform RO into PO. 

4 The inner circle 

The effective couplings € C, due to the imaginary parts of the diagrams. In the past this fact had 
some relevance only for realistic distributions while for PO they were conventionally taken to be real. 
The above convention has changed lately with the introduction of next-to-leading (NLO) corrections: 
imaginary parts, although not NLO in a strict sense, are sizeable two-loop effects. 

The explicit formulae for the Zff vertex are always written starting from a Born-like form of a 
pre-factor x fcrmionic current, where the Born parameters are promoted to effective, scale-dependent 
parameters, 



(if + i a L ) 7+ - 2 Q f 4 s 2 + w Q ]= 7m K + G{ 7s) , (6) 

where 7+ = 1 + 75 and a Q L are the SM imaginary parts. By definition, the total and partial widths of 
the Z boson include also QED and QCD corrections. The partial decay width is therefore described by 
the following expression: 

r f = r (z //) = 4 Cf r (\gf I 2 Rl + \g{\ 2 r{) + a ew/qcd , (?) 

where c/ = 1 or 3 for leptons or quarks (/ = I, q), and rC and i?{ describe the final state QED and QCD 
corrections and take into account the fcrmion mass mj. The last term, 

a -^r(') (8) 

EW/QCD x EW/QCD ^ EW ' \ U > 

accounts for the non-factorizable corrections. The standard partial width, r , is 

G M 3 

Tr " 82.945(7) MeV. (9) 



24\/2 7r 

The peak hadronic and leptonic cross-sections are defined by 



i27r S^ ^ 0=i27r ^> (io) 



z z z z 
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where V z is the total decay width of the Z boson, i.e, the sum of all partial decay widths. The effective 
electroweak mixing angles [effective sinuses) are always defined by 



4^^^=1-^4=1-4 



where we define 



The forward-backward asymmetry A is defined via 



(11) 
(12) 



.4. 



(13) 



where a F and er B are the cross sections for forward and backward scattering, respectively. Before analyzing 
the forward-backward asymmetries we have to describe the inclusion of imaginary parts. A FB is calculated 
as 

:i rr 

(14) 



A = 3 o " VA 

FB 4 cr^ 



where 



G F Ml 



V2 
8tt 



y/PePf QeQ/Re 



a*(M?)g e v g{x(s) 



PePfRe g% (px)* Re gf (g{)* s\x( S )\ 2 



(15) 



This result is valid in the realization where p/ is a real quantity, i.e., the imaginary parts are not re- 
summed in pf. In this case 



gf = Re(gl)+ibn(gt)=gf+iTm(g£), 

g{ = if + nm(g{). 



(16) 



Otherwise Q^ A — l^T 1 is a real quantity but pj is complex valued and Eq. (|l5|) has to be changed accordingly, 



i.e., we introduce 
with 



Pf v f 



(17) 



v f = lf -2Q fS m 2 9 f eS 



(18) 



For the peak asymmetry, the presence of p's is irrelevant since they will cancel in the ratio. We have 



FB A J J 



gin 



(19) 



The question is what to do with imaginary parts in Eq.(|19|). For partial widths, as they absorb all 
corrections, the convention is to use 



(R<J +(lm^ 



(20) 



On the contrary, the PO peak asymmetry A® f B will be defined by an analogy of equation Eq.(|l9|) where 
conventionally imaginary parts are not included 



A of — - A At 

A fb — A Steffi 



(g f v) 2 +(glY 



(21) 
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A definition of the PO heavy quark forward-backward asymmetry parameter which would include mass 
effects is 

2g b g b 

Ab = W^^m^^W 2 p ' (22) 

where (3 is the &-quark velocity. The difference is very small, due to an accidental cancellation of the 
mass corrections between the numerator and denominator of Eq.(p2|). This occurs for down quarks where 
(5t) 2 ~(9') 2 /2andwhere 



^ 3 2g%g\ 2g b v g\ 



m 2 q (9 b A ) 2 /2-(g b v r 

■s (gi) 2 + (g b A ) 2 



Therefore, our definition of the PO forward-backward asymmetry and coupling parameter will be as in 
Eq.©. 



5 Status of art 

The most important upgradings in the PO calculations consist of the inclusion of higher-order QCD 
corrections, mixed electroweak-QCD corrections || , NLO two-loop corrections of O (a 2 m 2 t ) @. The 
only case that is not covered is the one of final &-quarks, because it involves non-universal O (a 2 m 2 ^ 
vertex corrections. 

Another development in the computation of radiative corrections to the hadronic decay of the Z is 
contained in two papers Q , which together provide complete corrections of O (act s ) to T(Z — > qq) with 
q = u,d,s,c and b. 



6 Improved I/O Parameters 

Since all Renormalization Schemes (RS) use G F , we refer to the recent work of |t]] giving an improved 
value of G F = 1.16637(1) x 1(T 5 GeV~ 2 . An important issue concerns the evaluation of oqed at the 
mass of the Z. Define 

a(M z ) = ^ , (24) 

1 - Aa(5) (M z ) - A top (M z ) - A" Q P S (M z ) 

where one has 

Aq (5) [M z ) = Aa lept + Aa^. (25) 

The input parameter is now ActuA, as it is the contribution with the largest uncertainty, while the 
calculation of the top contributions to Aa is left for the code. Next, include for Aai op t the O (a 3 ) 

terms of § and use as default Ac^d = 0.0280398, taken from @. Therefore l/a {b \M z ) = 128.877, to 
which one must add the it contribution and the O (aa s ) correction induced by the it loop with gluon 
exchange [Tcj] . 



7 MI Calculus 

To summarize the MI ansatz, one starts with the SM, which introduces complex- valued couplings, cal- 
culated to some order in perturbation theory; next we define g' v , g^ A as the real parts of the effective 
couplings and T/ as the physical partial width absorbing all radiative corrections including the imaginary 
parts of the couplings and fermion mass effects. Furthermore, 

R q = ^, Ri = ^r, (26) 
1 h 1 1 
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for quarks and leptons, respectively. 

The experimental collaborations report PO for the following sets: 

(R f ,A^), (gf Y ,g{), (sin 2 <, P/ ). (27) 

In order to extract g^ , g^ from T/ one has to get the SM-remnant from Eq.(Q), all else is trivial. However, 
the parameter transformation cannot be completely MI, due to the residual SM dependence appearing 
inside Eq. ([?]). In conclusion, the flow of the calculation requested by the experimental Collaborations is: 

1. pick the Lagrangian parameters m t ,M H etc. for the explicit calculation of the residual SM- 
dependent part; 

2. perform the SM initialization of everything, such as imaginary parts etc. giving, among other 
things, the complement SM; 

3. select g f v ,g{; 

4. perform a SM-like calculation of T j, but using arbitrary values for g{,,g^, and only the rest, namely 



R f R f 



EW/QCD ' 



Im Ql , Im Q{ , 



(28) 



from the SM. 

An example of the parameter transformations is the following: starting from M Z ,T Z , R e „ T and 



j\0,e,n,T wc g rs ^. obtain 

FB 



With 



we subtract QED radiation, 



r e = m z t z 



and get 



12 TT R P 



1/2 



r u = 



r, = m z y z 



a R "1 1/2 



and 7 



12 7T 



6a/2tt ' 



1 + T 



3 a(M|) ' 

4 7T 



(29) 

(30) 
(31) 



sm 



I 1 



With 



we further obtain 



1 



2sin 2 # c e ff 



A f = - 



1 

4 + 
4 A°J 



3 »A e 



Im Q\ 



(32) 



(33) 



sin 6 



2 a f 



cir 



4 Qf 



-4/ 



r o r 



i -2|Q / |sin 2 ^ ff 



1 

4 



(im g{ 



-i -l 



(34) 



where / = u,t. For quarks one should remember to subtract first non-factorizable terms and then to 

r f 

distinguish between R v and R A . 



5 



Observable 


TOPAZO 


ZFITTER 


10 3 x T T Z 


l/a^(M z ) 


128.877 


128.877 




l/a(M z ) 


128.887 


128.887 




M w [GeV] 


80.3731 


80.3738 


-0.009 


< [nb] 


41.4761 


41.4777 


-0.04 


a" [nb] 


1.9995 


1.9997 


-0.12 


T h fGeVl 

lb 1 J 


1.74211 


1.74223 


-0.07 


[GeV] 


2.49549 


2.49573 


-0.10 


T„ [Me VI 


167.207 


167.234 


-0.16 


T e [MeV] 


83.983 


83.995 


-0.14 


T„ [Me VI 


83.983 


83.995 


-0.14 


T r [MeV] 


83.793 


83.805 


-0.14 


T„ MeV 


300.129 


300.154 


-0.08 


T d [MeV] 


382.961 


382.996 


-0.09 


T c [MeV] 


300.069 


300.092 


-0.08 


T h MeV 


375.997 


375.993 


0.01 


r inv [GeV] 


0.50162 


0.50170 


-0.16 




20.7435 


20.7420 


0.07 


#6 


0.215829 


0.215811 


0.08 


i?" 


0.172245 


0.172246 


-0.01 


• 2 zilcpt 

sm 0„ff 

en 


0.231596 


0.231601 


-0.02 


sin 2 $L 

en 


0.232864 


0.232950 


-0.37 


sin 2 d c cS 


0.231491 


0.231495 


-0.02 


Pe 


1.00513 


1.00528 


-0.15 


Pb 


0.99413 


0.99424 


-0.11 


Pc 


1.00582 


1.00598 


-0.16 


Observable 


TOPAZO 


ZFITTER 


10 3 x (T — Z) 


A li,l 

FB 


0.016084 


0.016074 


0.01 


A 0.b 

FB 


0.102594 


0.102617 


-0.02 


^B C 
FB 


0.073324 


0.073300 


0.02 




0.146440 


0.146396 


0.04 


A b 


0.934654 


0.934607 


0.05 




0.667609 


0.667595 


0.01 



Table 1: Complete table of PO, from TOPAZO and ZFITTER. 

8 Results for PO 

Having established a common input parameter set (IPS) we now turn to discussing the results for pseudo- 
observables (PO). The full list contains more PO and is given in Tab.(0), where we have included the 
relative and absolute difference between TOPAZO and ZFITTER in units of per-mill. 

9 Theretical Uncertainties for PO 

Here we discuss the theoretical uncertainties associated with PO. In Tabs. We give the central value, 
the minus error and the plus error as predicted by TOPAZO and compare with the current total experimental 
error if available. The procedure is straightforward: both codes have a preferred calculational setup 
and options to be varied, options having to do with the remaining theoretical uncertainties and the 
corresponding implementation of higher order terms. To give an example, we have now LO and NLO 
two-loop EW corrections but we are still missing the NNLO ones and this allows for variations in the 
final recipe for pf, etc. 
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TOPAZO has been run over all the remaining (after implementation of NLO) options and all the results 
for PO have been collected. We will use 

- central for PO evaluated at the preferred setup; 

- minus error for PO cenfra ; - min opt PO; 

- plus error for max opt PO - PO cenfra ;. 



O hsprva 1~>1 p 


Cfv\trf\] 


minim orror 


nlns prror 


total PYn frror 


L l a \ M z) 


1 OQ 077 
IZO.O / / 








1 In! M \ 


1ZO.OO / 








m w [(je\\ 


0U.0 / 01 


O.O IVlc V 


U.O 1V1C V 


yu ivic v 




41.4 / Dl 


1 n t^k 
i.u pu 


1.0 pD 


58 pb 


r„ [MeV] 


167.207 


0.017 


0.001 




T e [MeV] 


83.983 


0.010 


0.0005 


0.10* 


r M [MeV] 


83.983 


0.010 


0.0005 




T T [MeV] 


83.793 


0.010 


0.0005 




T„ [MeV] 


300.129 


0.047 


0.013 




T d [MeV] 


382.961 


0.054 


0.010 




T c [MeV] 


300.069 


0.047 


0.013 




T b [MeV] 


375.997 


0.208 


0.077 




Thad [GeV] 


1.74211 


0.26 MeV 


0.11 MeV 


2.3 MeV* 


Tinv [GeV] 


0.50162 


0.05 MeV 


0.002 MeV 


1.8 MeV* 


T z [GeV] 


2.49549 


0.34 MeV 


0.11 MeV 


2.4 MeV 



Table 2: Theoretical uncertainties for PO from TOPAZO. *) assumes lepton universality. 



10 Realistic Observables 

The RO are computed in the context of the SM, however one of the goals will be to pin down the definition 
of PO; the calculation of RO in terms of the defined PO for the purpose of MI fits, showing that for PO 
with values as calculated in the SM, the RO are by construction identical to the full SM RO calculation. 
The last point requires expressing p's and effective mixing angles in terms of PO, assuming the validity of 
the SM. One should remember that gauge invariance (GI) at the Z pole (on-shell GI) is entirely another 
story from GI at any arbitrary scale (off-shell GI). 

Some of the re-summations that are allowed at the pole and that heavily influence the definition of 
effective Z couplings arc not trivially extendible to the off-shell case. Therefore, the expression for RO 
= RO(PO), at arbitrary s, requires a careful examination and should be better understood as 

RO = RO (PO,SM) , (35) 

that is, for example: 

cmi = cr SM (Ri,A°' B , 



g f v ,g{ 



Pf, sin 2 6^ s ; residual SM^j . 

(36) 
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Observable 


central 


minus error 


plus error 


total exp. error 


Ri 

R° c 


20.7435 
0.215829 
0.172245 


0.0020 
0.000100 
0.000005 


0.0013 
0.000031 
0.000024 


0.026 
0.00074 
0.0044 


sin 2 ^ 
sin 2 e b ff 
sin 2 ^ ff 


0.231596 
0.232864 
0.231491 


0.000035 
0.000002 
0.000029 


0.000033 
0.000048 
0.000033 


0.00018 


A o,i 

4i 

A 0,c 

FR 


0.016084 
0.102594 
0.073324 


0.000057 
0.000184 
0.000142 


0.000060 
0.000195 
0.000149 


0.00096 
0.0021 
0.0044 


A h 
A c 


0.146440 
0.934654 
0.667609 


0.000259 
0.000032 
0.000114 


0.000275 
0.000005 
0.000103 


0.0051 
0.035 
0.040 


Pe 
Pb 
Pc 


1.00513 
0.99413 
1.00582 


0.00010 
0.00048 
0.00010 


0.000005 
0.000001 
0.000005 


0.0012 



Table 3: Theoretical uncertainties for PO from T0PAZ0. 



As long as the procedure does not violate GI and the PO are given SM values, there is nothing wrong 
with the calculations. It is clear that in this case the SM RO coincide with the MI RO. 



11 Final-State Radiation 

One should realize that s' (s' being the centre-of-mass energy of the e + e~ system after initial state 
radiation) is not equivalent to the invariant mass of the final-state // system, M 2 (ff), due to final state 
QED and QCD radiation. Furthermore, in the presence of a s'-cut the correction for final state QED 
radiation is simply 

r fs _ 3 n2 <*(*) 
QED - 4 Hf , {o<) 

while for a M 2 -cut the correction is more complicated, see pT| . 
For full angular acceptance one derives the following corrections: 



4-7T 



— + 2 In x ) In -4y 

2 J mj 



Here we have introduced 



-z(l + -) Inz + 2C(2) - 2Li 2 (x) - 2hxx + - - 3z - - 
2/ 4 ' 

x=l-z, z = M 2 (//) /s. (38) 

For an s'-cut both QED and QCD final-state radiation are included through an inclusive correction factor. 
For M 2 -cut the result remains perfectly defined for leptons, however for hadronic final states there is a 
problem. This has to do with QCD final-state corrections. Indeed wc face the following situation: 

1. for e + e~ — * ffj the exact correction factor is known at O (a) even in the presence of a M 2 (ff) 
cut; 
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2. the complete set of final-state QCD corrections are known up to O (a 3 ) |T^] only for the fully 
inclusive setup, i.e., no cut on the // invariant mass; 

3. the mixed two-loop QED/QCD final-state corrections are also known only for a fully inclusive 
setup Q 

4. at O (a s ) QCD final-state corrections in presence of a M 2 -cut follow from the analogous QED 
calculation of . 



12 Convoluted Realistic Observables 

We have the following options: 

CA3 Complete RO, with QED initial state radiation implemented through an additive O (a 3 ) radia- 
tor || 

CF3 Complete RO, with QED initial state radiation implemented through a factorized O (a 3 ) radia- 
tor @ 

The following two equations define cross-sections and forward-backward asymmetry convoluted with 
ISR: 

ct t ( s ) = / dzH (z; s) <t t (zs), (39) 

where z a — sq/s and 

Tra 2 QlQ 2 f f 1 1 

4b ( s ) = / dz ~, — 75 H FB ( z ; s ) ^fb (**) , (40) 

°tot Jz (1 + Z) 

Note that the so-called radiator (or flux function), H(z;s), is known up to terms of order a 3 while 
H FB is only known up to terms of order a 2 . The kernel cross-sections tr T FB should be understood as the 
improved Born approximation (IBA), corrected with all electroweak and possibly all FSR (QED (g> QCD) 
corrections, where all coupling constants and effective vector and axial weak couplings are running, i.e., 
they depend on s' — zs. 



13 Uncertainty on QED Convolution 

We define the following quantities: 

SdeC (°) = 0$5 ~ 1 > A dcc (0) =0-O su . (41) 

For convenience of the reader we have reproduced in Tab.(||) the results for the CA3 and CF3 mode. 
From Tab.(Q) we derive the absolute differences, for A^ B , and the relative ones, for cross-sections, between 
additive and factorized versions of the QED radiators. They are shown in Tab.(|h. 



14 Initial-Final QED Interference 

If initial-final QED (ISR-FSR) interference (IFI) is included in a calculation we have a conceptual problem 
with the meaning of the s'-cut. In this case the definition of the variable s' is unnatural since one does 
no longer know the origin of the radiative photon (ISR or FSR). Only a cut on the invariant mass of 
the final-state // system would make sense. There is another option: to select events with little initial- 
state radiation one can use a cut on the acollinearity angle 9& C oi, between the outgoing fermion and 
anti-fcrmion. A cut on # a coi is roughly equivalent to a cut on the invariant mass of the // system, indeed 
one may write 

_ 1 - sin(fl acol /2) r49 s 
s eff ~ l + sin(W2) ' ( ' 
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therefore a cut of a coi < 10° is roughly corresponding to the request that s'/s > 0.84. The inclusion of 
initial- final QED interference in TDPAZO and in ZFITTER is done at O (a), i.e., the O (a) interference term 
is added linearly to the cross-sections without entering the convolution with ISR and without cross-talk 
to FSR. For loose cuts the induced uncertainty is rather small and the effect of the interference itself 
is minute. On the contrary, when we select a tight acollinearity cut the resulting limit on the energy 
of the emitted photon becomes more stringent and the effect of interference grows. The corresponding 
theoretical uncertainty, due to missing higher-order corrections, is therefore expected to be larger. 

I/F INT is currently treated dubiously by the experiments: The MC used to extrapolate for efficiency 
and possibly acceptance do not contain I/F interference (eg, KORALZ in multi-photon mode). Thus 
extrapolated and quoted results somehow miss I/F interference. There is a claim of some serious discrep- 
ancy between KK [|l7) / T0PAZ0 / ZFITTER in cr int , for muons (see frittp: / /home. corn. ch/jadacb| ) It has 
to do with exponentiation of ISR (8 FSR. 





LEP 1 energy in GeV 


M z -3 


M z - 1.8 


M z 


M z + 1.8 


M z +3 


Jdecjo-CAS) z 
£dec ((7 CF3) T 

6 dcc (a^ F3 ) Z 


-23.976 
-24.007 
-23.973 
-24.000 


-27.616 
-27.629 
-27.611 
-27.624 


-26.257 
-26.261 
-26.253 
-26.256 


5.356 
5.341 
5.364 
5.348 


30.665 
30.631 
30.671 
30.637 


<5 deC (^a A d 3 ) Z 

<5 d ° C K C a F d 3 ) T 
^Wd 3 ) Z 


-25.867 
-25.873 
-25.863 
-25.869 


-28.501 
-28.503 
-28.497 
-28.499 


-26.537 
-26.538 
-26.533 
-26.533 


4.952 
4.945 
4.959 
4.953 


30.564 
30.550 
30.572 
30.558 


A dcc (Ap A3 ) Z 

^dec^/iCF3^ j 
A dcc (A^F3) z 


-2.190 
-2.211 
-2.189 
-2.215 


-1.967 
-1.975 
-1.966 
-1.977 


-1.804 
-1.803 
-1.804 
-1.803 


-6.292 
-6.287 
-6.292 
-6.286 


-11.486 
-11.480 
-11.487 
-11.479 



Table 4: RO: the effect in % of initial state QED radiation for CA3 and CF3 modes. 





LEP 1 energy in GeV 


M z -3 


M z - 1 M z M z + 1 


M z +3 


10 4 x (fact/add-1) 


C M 


0.44 
0.88 


0.63 
0.63 


0.61 
0.68 


0.72 
0.72 


0.49 
0.49 


Chad 


0.58 
0.61 


0.58 
0.62 


0.64 
0.67 


0.73 
0.76 


0.59 
0.62 


fact- add [pb] 


c M 


0.01 
0.02 


0.03 
0.03 


0.09 
0.10 


0.05 
0.05 


0.02 
0.02 


Chad 


0.26 
0.27 


0.56 
0.60 


1.95 
2.04 


1.04 
1.08 


0.48 
0.51 


10 5 x (fact-add) 


A£ 

FB 


1.00 
-4.00 


1.00 
-2.00 


0.00 
0.00 


0.00 
1.00 


-1.00 
1.00 



Table 5: Absolute and relative differences in TDPAZO and in ZFITTER for additive and factorized radiators. 
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15 Note added: upgrading of ISPP 



Here we discuss some recent work connected with Initial State Pair Production (ISPP). We consider two 
different approaches: the first is based on JMS [l8| |, an ansatz for re-summation of both soft photons and 
soft lepton pairs. For hadron pairs it is corrected numerically with an uncertainty of 30%. 

All versions of TDPAZO < 4.4, use an effective, naive, implementation of ISPP (ONP = Y) based 
on KKKS jl9) interfaced with soft photon re-summation. The agreement with JMS is perfect around the 
resonance but quickly deteriorates on the high energy side. 

In the latest version of TDPAZO, the variable ONP has the additional value of ONP = I, where 

• KF approach p0| for virtual pairs and for soft and exponentiated pairs is used; it is also applicable 
to hadron pairs where we use KKKS results for O (a 2 ) and write them in terms of moments. Then 
we match it to KF and generalize KF to soft and exponentiated hadrons pairs. 

• Next, we use generalized KF for virtual and soft pairs and cut to the same s' cut of IS QED 
radiation. Once an s' cut is introduced in one place then it should be used everywhere, even for 
photons + pairs. 

• Finally, we use the soft approximation only up to some cut A that is compatible with £> A>2m. 
Above the cut we use KKKS-formalism but not added linearly to the cross-section; instead we use 
a KKKS correction factor in convolution with IS QED radiation. The radiator used here is a 
leading- log (LL) radiator, evaluated at the second order or, optionally, at (third order). 

Now, some comparison^, the cross-sections at the Mi-point M z = 91.1882 GeV, T z = 2.4952 GeV, 
a\ = 41.560 nb and R x = 20.728, A™ = 0.017319, for ISPP0 on and off; ISPP effect was taken as the 
difference. We report Chad [nb] in Tab.(|^): first entry is without pairs, second entry is with pairs and 
third entry is (with - without). 

Next, some SM numbers, i.e. the hadronic cross-section for M z = 91.1867 GeV, mt — 173.8 GeV, 
M H = 100 GeV and a s (M z ) = 0.119. In Tab.(0) first entry is ONP = N, second entry is ONP = I and 
third entry is with/without - 1 in per-mill. 
Similar results have been obtained in fl21fl . 

For an extension to LEP 2 energies [E2[ let's start from a simple case, e + e~ PP-corrections to e + e~ — > 
66. To avoid confusion we will introduce the following definition: Multi-Peripheral or MP diagrams; 
Initial State Singlet, or ISS diagrams; Initial State Non-Singlet, or ISNS diagrams; Final State, 
or FS diagrams. 

Note that we include both 7 and Z exchange, so that one could still distinguish between ISNS 7 , 
ISNSz and interference. On top of real pair production one has to include virtual e + e~ pairs. 

So far, TOPAZO (and ZFITTER too) only include ISNS 7 plus virtual, with ISS optional. The approxima- 
tion is well justified around the Z-peak but now we have to move to higher energies, where new thresholds 
open. Furthermore, one has to add a proper definition of soft pairs and of hard pairs. We can define 
Soft(Hard) Invariant Mass pairs, or SIM(HIM) pairs, according to some pair mass cut. Of course, one 
could still use a soft-hard separation based on other variables but, in any case, also very hard SIM pairs 
are 2f signal. 

We need an operative, universally accepted, separation of various contributions. Consider again the 
process e + e~ — > 66e + e _ ; the total cross-section contributes to three different processes: a) genuine 4f 
events; b) e + e~ PP-correction to e + e~ — > 66; c) 66 PP-correction to Bhabha scattering. 

A naive separation, easy-to-implement, would be the following. Let the process be specified by 
e+e~ — > bb(Q 2 ) + e + e~(q 2 ). We can integrate and reduce the cross-section to a two-fold integral with 
the following boundaries: 

4m 2 < Q 2 < (V« - 2m e ) 2 , 

4 m 2 < q 2 < VO 2 ) 2 . (43) 

Next, we introduce two cuts z p , z s , i.e. primary and secondary cuts and define 

2 JMS numbers are a courtesy of B. Pietrzyk and of G. Quast 

3 Here ISPP means NON-SINGLET Initial State Pair Production. SINGLET pairs are, in principle, subtracted 
together with two-photon contributions 



11 



Vi[GeV] 


JMS 


T0PAZ0(ONP=Y) 


T0PAZ0(ONP=I) 


88.464 


5.2297 


5.2187 


- 




5.2170 


5.2059 


5.2063 




-0.0127 


-0.0128 


-0.0124 


89.455 


10.121 


10.105 


- 




10.094 


10.079 


10.079 




-0.027 


-0.026 


-0.026 


90.212 


18.238 


18.227 


- 




18.187 


18.178 


18.179 




-0.051 


-0.049 


-0.048 


91.207 


30.529 


30.549 


- 




30.449 


30.471 


30.470 




-0.080 


-0.078 


-0.079 


91.238 


30.589 


30.609 


- 




30.510 


30.531 


30.530 




-0.079 


-0.078 


-0.079 


91.952 


25.173 


25.166 


- 




25.121 


25.117 


25.112 




-0.052 


-0.049 


-0.054 


92.952 


14.503 


14.484 


- 




14.488 


14.476 


14.466 




-0.015 


-0.008 


-0.018 


93.701 


10.067 


10.051 






10.064 


10.059 


10.046 




-0.003 


0.008 


-0.005 



Table 6: Comparison of T0PAZ0 with JMS-approach for the inclusion of ISPP. 



V^[GcV] 


s'/s = 0.01 


s'/s = 0.1 


s'/s = 0.5 


s'/s = 0.9 


M z - 3 [GeV] 


4.45106 


4.44444 


4.43769 


4.38335 




4.44069 


4.43391 


4.42705 


4.37227 




-2.33 


-2.37 


-2.40 


-2.52 


M z - 1.8 [GeV] 


9.60037 


9.59393 


9.58742 


9.52026 




9.57589 


9.56929 


9.56267 


9.49497 




-2.55 


-2.57 


-2.58 


-2.66 


M z [GeV] 


30.43633 


30.43005 


30.42316 


30.31927 




30.35752 


30.35109 


30.34409 


30.23944 




-2.59 


-2.59 


-2.60 


-2.63 


M z + 1.8 [GeV] 


14.18153 


14.17541 


14.16813 


13.98373 




14.16470 


14.15846 


14.15117 


13.96556 




-1.19 


-1.20 


-1.20 


-1.30 


M z + 3 [GeV] 


8.19803 


8.19206 


8.18498 


7.87534 




8.19800 


8.19193 


8.18460 


7.87272 




-0.004 


-0.02 


-0.05 


-0.33 



Table 7: SM results with the inclusion of ISPP. 
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• e + e PP-correction to e + e — ► bb: 

z p s < Q 2 < (y/s - 2m e ) 2 , 4ml <q 2 <min{z s s, (y/s - \/0 2 ~) 2 }. ( 44 ) 

• bb PP-correction to Bhabha scattering: 

z p s < q 2 < (y/s- 2m b ) 2 , 4m 2 b < Q 2 < min{z s s , (y/s - x/q 2 ) 2 }. (45) 

• while the remaining portion of the phase space is background. 

At the same time, we must address what has to go in the calculation and what is already subtracted 
from the data. We have to define the exact content of PP-corrections to be inserted in the upgradings 
of the semi-analytical codes. Then, one can apply a theory correction, actually a subtraction obtained 
from MC, to go to signal definitions as employed by TOPAZO/ZFITTER. This is what is done now, the only 
problem being that the present signal definition in TOPAZO/ZFITTER is not good enough for LEP 2. Thus 
signal definition must be improved and we have to find an agreement on the optimal signal/background 
separation; it should be realistic enough to be used in the experimental analysis but simple enough to 
allow its implementation in the codes. 

Typical example? Using invariant-mass cuts one finds relatively large corrections (around 5% or 
more), most of which ares due to the MP process. So, if the situation remains unaltered, one has to apply 
a rather large subtraction, which is rather unnatural. Of course, we can take the MP into account, but 
only if some reasonable cut is applied to M(bb). At 200 GeV the double resonant e + e~ — ► ZZ — > e + e~bb 
cross-section, i.e. the NC02 MC, is about 10 _2 pb, while the total multi-peripheral, that is not accessible 
with a NC02 (two diagrams) MC, is about 1 pb compared with a <7h a d of 80 — 90 pb. 

However, this huge increase in cross-section is caused by hadronic two-photon collision processes, 
where the bb system is created by the two photons radiated of the e + e~. Thus the invariant mass of the 
bb system will be very low, peaking toward 2mj. The bulk of the cross section will be around 10 GeV 
invariant bb invariant mass. 

Finally we come to the most complicate configurations, e + e~ — > qqQQ or even worst, e+e~ — > qqqq, 
e.g. uuuu etc. configurations. We have to generalize our previous separation of the 2f signal from the 
4f background. There are two options: we require that at least one pair has an invariant mass greater 
than z p s. But only the qq pairs or, should we also include the qQ pairs? And what about the qQ,qQ 
pairs? Or, we require that at least one pair has an invariant mass greater than z p s while all remaining 
configurations have invariant mass less than z s s. 

Now we come to the really difficult part of the problem. Since arbitrarily low SIM pairs are allowed, 
in both cases, a parton-level calculation cannot be accurate enough. Under the simplified assumption 
that one pair (the Q 2 one) is HIM enough and that the remaining one (the q 2 one) is SIM enough we 
can write the cross-section as 

do-f fa\ 2 , n2 .,Rh a d(q 2 ) 



U "/W)-3^3-/(«>0'O 



alq 2 dQ 2 \ttJ ' 3sq 

ft rfl 2\ s + (Q 2 + g 2 ) s-Q 2 -g 2 + x/\ r- 

/(S ' Q ' 9) = s-Q 2 -q 2 lD s -Q 2 - g2 -VA ' 2 ^ 

A = X(s,Q 2 iq 2 ). (46) 

Clearly we face two problems, one conceptual and another practical. The first problem will be referred to 
as the double- counting problem. At the parton level, we have e + e~ — > qqQQ and, for instance we want to 
sum over q in order to define the Q-line shape. However, when we also sum over Q in order to have the 
full hadronic line shape, proper care must be taken in order to avoid double-counting. At the parton level 
this can be done but one should remember that we actually compute e + e~ — ► QQ + hadrons, through 
-Rhad(s), and the problem is, by far, more severe. Therefore, we have to agree on what should happen 
once -ffhad(s) is called for low s. 

At least in principle, we can introduce three different cuts, i.e. s'/s for initial state QED radiation 
and z p ,z s for pair production. From a theoretical point of view, there is no problem here. We should 
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reach an agreement: one common cut (s'/s = z m [ n ), two cuts (s'/s z min but also s'/s = z p and z s ) or 
three cuts (s'/s z p z s)- 

Tentative conclusions are as follows: the whole 4f must be included to compute the 2f cross-sections; 
Or, the whole 4f is to be divided into two components, signal and background. For our purposes their 
definition is peculiar, signal is what you have implemented into the semi-analytical codes. Background is 
what one subtracts by using a MC. We go from one extreme solution to the other: 

Background — if everything is included in the semi-analytical calculation. Multi-peripheral is an 
example of something difficult to implement into the semi-analytical approach if low-invariant mass 
regions are required. 

Signal = ISNS, i.e. everything else (large effects) is subtracted by MC. However, using different MC 
programs would bring to subtractions that differ by some per-cent, which then would have to be regarded 
as a theoretical systematic uncertainty. 

One could decide to count as background 4f processes from true e + e~ annihilation (e.g. ZZ, WW) and 
4f processes via e-gamma collisions (e.g. (e)eZ, (e)vW). This distinction is based on selecting sub-sets of 
diagrams to be used for the 2f signal definition. However, we should not define the signal according to the 
kind diagrams it originates from. We need a definition based on measurable quantities, so we should not 
exclude a set of diagrams, we should exclude the kinematical region where that particular set dominates. 
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